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Abstract 

Direct analytical and numerical calculation show that two-electron atomic config- 
uration can be unstable with respect to a static or dynamic shift of the electron 
shells. This enables to develop a so called non-rigid shell model for a partial account 
of the electron correlations within atomic clusters in solids. In a framework of this 
model a correlated state of two-electron molecular configuration is described by a 
set of symmetrized shell shifts q'^ similarly to the well known shell model developed 
for a description of the lattice dynamics. A set of g^-shifts are found after mini- 
mization of the energy functional. We present a number of the novel unconventional 
effects including: i) a correlational mechanism of the local pairing; ii) a correlational 
(pseudo) Jahn- Teller effect provided by a joint account of the electron shell shifts 
and conventional nuclear displacements; iii) an appearance of the chiral correla- 
tional states. The model allows an introduction of the pseudo-spin formalism and 
effective "spin-Hamiltonian" for a description of the short- and long-range ordering 
of non-rigid atomic backgrounds in crystals. Finally, the model can be readily built 
in the conventional band schemes. 



Introduction 



Electronic correlations is one of the fundamental problem in a theory of atoms, 
molecules and solid state, particularly for the systems with high density of ex- 
cited states when a small perturbation can result both in drastic reconstruction 
of the energy spectrum and in modification of the ground state up to forma- 
tion of a strongly correlated state. As a rule, in such a situation an appropriate 
description of the ground state within the bare restricted basis often requires 
a lot of configurations or considerable extention of the basis, and so becomes 
difficult for practical realization and interpretation. Namely this situation oc- 
curs in atoms, where description of some specific correlation effects in terms 
of Hartree-Fock basis requires a large number of Hartree-Fock configurations. 
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Such a problem implies a search for alternative variational approaches to the 
electronic structure and energy spectrum. 

In this work we develop further the variational method for the many-electron 
atomic clusters with the trial parameters being the coordinates of the center 
of the one-particle atomic orbital [1]. The resulting shift of the atomic orbital 
allows to interpret the variation of the electronic density distribution rather 
clearly, and the symmetry of a system can be readily used for construction of 
the trial many-electron wave function. The shifted electronic shells in conven- 
tional MO-LCAO-scheme with restricted set of the one-particle states allow to 
take into account an additional multipole interactions and to construct novel 
states with unique properties. As a whole the model bears a strong resem- 
blance to the well known shell model widely used in lattice dynamics. 



2 Two-electron configuration 



We consider the problem of two electrons in certain atomic potential to be a 
simplest model for manifestation of electronic correlations. The orbital part 
of the singlet two-electron wave function formed by the shifted one-particle 
orbitals (bi-orbital) can be written as follows: 



fi,f2]a,(3) = T] ^ ^p{fi-a)i){f2-(^)+i){fi-i5)-i}){f2-a) , (1) 



where a, j3 are the displacement vectors for the one-particle orbitals (Fig.l), 
7] the normalization factor. Below, only the real functions of s-type are used 
as the trial one-particle states. Then 

77^ = 2 (1 + ^2 (a, (2) 



where 6'(a, /3) is the overlap integral for the one-particle orbitals. The Hamil- 
tonian of the problem in atomic units (eo = ^ = — , oq = -^) is 
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The variational procedure is performed with the full energy functional: 



= I ^ I = £; (a, p) . 



(4) 



Taking into account the expression (1) we obtain: 
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E{d, P) = [2 d) - Zo a) + u0, + (5) 

2 /3) /3) - 2Zo S{d, /3) /3) + c(a, /3) + a(a, 

where the following matrix elements are introduced: the one-center integral 
a) is the kinetic energy of an electron with functions of the same cen- 
ter, the two-center integral a) is the interaction of an electron with the 
potential center with functions of the same center, the two-center integral 
(3) is the kinetic energy of an electron with functions of different centers, 
the three-center integral u{d, /3) is the interaction of an electron with the po- 
tential center with functions of different centers, c(a, /3) and a(a, /3) are the 
Coulomb and the exchange parts of the inter-electron interaction. One should 

— * — * — * 

note that S{d,P), c{d,P), a{d, P) are the two-center integrals. 

In further we use the Slater functions with the index k and effective charge Z 
as the one-particle atomic orbitals: 

V'(r)-^^z,ferV^^ (6) 
where the normalization factor is 



The analytic expressions for the matrix elements are presented in Appendix. 
Their examination allows to obtain some information about extremal values 
of the d and p. Introducing 

q+ = ^(d + P),q. = ^(d-p) (8) 

and the coordinate system with center at g+ = 0, one can see that only 
u{d, d) + u{P, /?), u{d, (3) depend on 

u{d,d)+u0J) = I (^2(f-g_)+t^2(f+g_)), (9) 

u(d, /3) = / I ^ I i/jif— g_) q-). 
J I r — g+ I 

These quantities are invariants with respect to inversion in the displacement 
vector space, hence = is a critical point in g!|_-space. The surface E{'^} = 
const in g_|_-space is a sphere for the one-particle s-functions , and the point 
= is a minimum or a maximum. In general, for an arbitrary angular 
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dependence of the one-particle functions ijj the point q+ — can be also a 
saddle point. 

It is rather clear, why the point q+ — appears to be critical: at the given value 
of the interaction with the potential center, the configuration with a = —j3 
can minimize inter-electron repulsion. This is confirmed by results of numerical 
minimization of the full energy functional for the Is^-configuration {k = 0) 
which are listed in the Table 1 and 2. Thus, the following function can be 
defined 

f2; q) = 77"^ [ip{fi - q) ip{f2 + q) + ip{fi + q) ip{f2 - q)] , (10) 

that leads to the reduction of the number of trial parameters. The full energy 
functional with functions ip{r) (6) depends only on the q = \q\ , but not on 
the direction of q. 

The value Z in the calculations mentioned above was a free parameter. This 
parameter provides an additional mechanism of the electronic density redistri- 
bution along with the electronic shell shift. For the isolated atom it is naturally 
to assume that Z is also the variational parameter, because in this situation 
the atomic potential is the only mechanism of the electronic density redistri- 
bution: 

E m = E{q, Z) = E{a ^ qj ^ -q; Z) (11) 

It should be noted that in a crystal the functions with certain " atomic" value 
Z which minimizes the energy of the isolated atom can form the strongly 
correlated state like shifted electronic shell state for the minimization of given 
crystal potential characterized by the parameter Zq. 



3 Expansion of the full energy functional 



The function (10) possess the following property: it has no linear term in the 
g-expansion at g = 0. As 



dip{f— q) d%l){f 



q=0 — ^ |g=0) 



dq dq 



(12) 



the first derivative of the function (10) at g = turns to zero. Hence, the 
E{q, Z) with the functions (10) has an extremum at g = 0, which type is 
defined by sign of the quadratic term E^"^^ in the g-expansion of the -E'(g): 

E(g) g". (13) 
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For the functions (6) the full energy functional doesn't depend on the direc- 
tion of (f, and therefore we consider shifts to be directed along ^-axis. The 
g-expansion of the function (10) can be written in the following form: 

■0(f — q) ~ a(r) — 6(r) q H — c(r) g^, (14) 



where 



a{r) = V'(r - q) \,=o = Nz,k r'e-^^ (15) 
b{r-)^^^p^U=o^Nz,k{kr'-'-Zr') e-^^cos^, (16) 

^(^1 - ^ \i=o = NzAkr'-' - Zr'-' + (17) 

+ {kik - 2) r*^-' - Z{2k - 1) r''-' + Z'^r') cos' 6] e-^^ 

with the angle 9 counted out from 2;-axis. The g-expansion of the matrix 
elements up to quadratic terms and the expressions for E'^^^ and E^"^^ are 
listed in the Table 3. 

From the expression for E'^'^'> the criterion of the non-zero shift of the electronic 
shell can be obtained for the Is-function {k = 0) as the one-particle state. The 
shift is not zero, if Z > — in opposite case there is no displacement of 
the electronic shell from the potential center. This result is entirely compatible 
with our numerical calculations presented in Table 2. The minimum of E{Z^ q) 
for /c = is obtained for g^m = and Zmin = — jq, that agrees with the 
well-known result in the helium atom theory [3] . 

From the other hand, the full energy functional with functions (10) formed 
from the one-particle ns-states (6) with k ^ has the minimum at g 7^ for 
any values Zq and Z . It means that for these states the non-zero shift of the 
electronic shell always takes place. 

The expression for E^^^ allows to find the value of the effective charge Z that 
minimizes full energy of the ns'-configuration at g = 0: 

_^2fc + l ( 2k + 1 {2k + l)(Ak + 3)\ \ 
°)t+l [2k + 2 24fc+3[(2A; + 2)!]V' 



At A; ^ 00 this expression tends to Z^^^ — 2Zo — 1. 

Origin of the different behaviour of the electronic shells with k = and k 7^ 
can be readily understood from the data listed in Table 3. The principal 
difference of the one-particle function with k ^ from that of with k — is 
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that the function with k ^ turns into zero at r = (Fig.2). The energy of 
interaction of electrons with the potential center appears to be most sensitive 
to the electron density on nucleus. Its displacement leads to strong increase 
of the full energy of the system in the case of wave function with k — 0, and, 
from the other hand, zero density at r = for the function with k ^ provides 
specific "softness" of this part of interaction that is revealed in zero value of 
quadratic term in expansion of u{q, —q)- The expressions for matrix element 
^(^i ~Q) describing relative motion of the electron shells are different at A; = 
and A; 7^ whereas the overlap integrals («S'(g, — g), c(g, — g), a(g, —q)) have the 
same structure for the both cases. Finally, it can be concluded that the gain in 
the electron-electron interaction at A; 7^ with displacement of the electronic 
shells hasn't been compensated by the loss in the energy of interaction with 
the potential center as in the case A; = 0. 

An appearance of the electronic density on the potential center also implies the 
gain in energy. This makes possible to determine the most favorable directions 
of displacement of the electronic shells in the case of the anisotropic wave 
function with node at r = 0. So, for p^-orbital it has to be z-axis direction, 
for (i^2_j^2-orbital the critical directions are along x- and |/-axis, and so on. 
Finally, for the A; 7^ orbitals one should expect a well developed non-trivial 
q ^ minimum at the energy surface E{q,Z). This conclusion appears to 
be compatible with the results of numerical minimization of the full energy 
functional for np^-configuration [4]. 



4 The form of the full energy functional E{q, Z) 

The general expressions for the matrix elements in case of the electronic con- 
figuration with a = q and P = —q can be obtained from those listed in Table 
1 at p = '^Zq, a — Zq and with transition to the limit ^ — > 0, 77 — > in 
expression for u{a, (3). In terms of p = 2Zq the latter can be written as: 

2fc+2 fc+1 

""^^^ ^ ^W^. + ^%^"\P)- (19) 

At A; = the matrix elements coincide with the well-known fundamental 
results from the atomic theory [5,6]. 

The results of minimization of E{q, Z) for a number of the lowest values of 
k {k = 0, ...4) are listed in Table 4; the full energy functional as a function 
of q at Z = is shown in Fig. 3. As it was mentioned above, for A; 7^ 

the displacement of the electronic shells takes place. The gain in energy is 
0.22-=- 0.36 a.ii., and the value of displacement is 0.34-=- 0.76 a.w. for k — 1, ..A. 



6 



The /c-dependencc of the global minimum can be explained by the character 
of the electronic density distribution at different A;. The Fig. 2 shows that the 
increasing in k is accompanied by the lower values of the electronic density 
near r — 0. The gain in energy of a system can be provided by decreasing of 
electron-electron interaction (decreasing of positive contribution from t{q, —q), 
O'iQT—Q)) and by increasing of the interaction with the potential 
center (increasing of negative contribution from u{q,q), u{q, —q))- The rising 
of qmin with the index k is associated with the delocalization '^'^(r); the lower 
value of the overlap and, therefore, the lower interaction of the electrons are 
obtained with rising of q. Another mechanism providing the overlap decrease 
is in increasing of Z, that leads also to a localization of il){r) (see Fig. 4). With 
increasing Z the negative contribution from interaction with the potential 
center and the positive ones from kinetic energy are increased simultaneously 
resulting in certain compromise value of Z. This value increases with k due 
to a delocalization of ■0(r) with k. In Fig. 5 the value qmin as function of Z is 
also shown. 



5 The dynamic shifts of the electronic shells 

The distribution of the electronic density at A; = 1 with and without shifts is 
shown in Fig. 6; in the both cases the wave functions providing the minimum 
of the full energy functional are used. The symmetry of the electronic density 
distribution with the shifts (Coo/i) breaks the initial spherical symmetry, which 
can be restored with taking into account the energy equivalent configurations 
with various directions of the displacement vectors. 

The full energy functional can have the continuum of the equivalent minima 
in the displacement vector space. In this sense the system has the variational 
degeneracy. The existence itself, the form and other parameters of the minima 
continuum depend on the one-particle states and the parameters of the po- 
tential. In the case of the ns^-configuration mentioned above, only the value 
of difference of the displacement vectors of the one-particle states is fixed 
\p\ = \q. — I3\ = 2q (with a — q, 13 = —q), so the continuum is a sphere in the 
q space that restores the initial spherical symmetry of the system. By analogy 
with the description of collective motion in nuclei [7] , one can form the linear 
combinations with a help of the bi-orbitals ^ which have different vectors qmin'- 

(ri, fa) = I ^ (fi, fa; q, -q) f (n) dn, (20) 

where integration is performed on a sphere in the q space. Such a linear com- 
bination can provide the lower energy due to the "off-diagonal" in q matrix 
elements of the full energy functional which take account of the "interaction" 
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of bi-orbitals. The variational procedure with the functions (20) yields an in- 
tegral equation for the function / {fl) : 



I dnf (fi) [K {q, q')-E-I (g, g')] = 0, 



(21) 



where 



K{q, cf) = {^{q, -q) \ H \ m{<f, -cf)) 
I{q,q') = {^{q,-mi^,-^))- 



With account of the symmetry of the ns^-configuration [7] the following solu- 
tions of (21) can be written: 



In other words, for the ns^-configuration with shifted shells the set of orthog- 
onal states can be introduced: 



the terms of which transform according to irreducible representations of the 
rotation group. These states could be called the dynamic ones, as they can 
result in the correlational contribution to orbital current. The spectrum of 
these states can be not similar to that of the space rotator. Note that for 
such states with the dynamic shifts of the electronic shells one might expect 
the anomalously large values of the electric (dipole, quadrupole) or magnetic 
susceptibilities, and these values reflect the electronic correlation effects. 



6 The MO-LCAO-scheme with the shifted atomic orbitals 

Introduction of the shifted one-electron atomic orbitals implies the general- 
ization of the conventional MO-LCAO-scheme [8]. Instead of standart set of 
the molecular orbitals (MO) iproyQ{f, 0) being the symmetrized combination of 
the atomic functions centered in the points of the equilibrium nuclei positions 
{qr-y — 0), the new set of the shifted MO: 



should be introduced, where qr-y is the symmetrized coordinate of the atomic 
orbital displacements in cluster, Tg^^ the operator of the symmetrized displace- 
ment. Such an approach is the natural generalization of the shifted electronic 



f{n)^YLM{o,v). 



(22) 




(23) 



(24) 
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shells model for the many-atomic cluster. The symmetry group of the wave 
function (24) is an intersection of the group germs of the Tq and G. In contrast 
with the symmetrized coordinates of the cluster vibrations, the vector q-p^ is 
fixed and defines the certain distorted distribution of the electronic density. 
Ur Ai, then the function (24) doesn't possess the proper transformation 
properties or, in other words, doesn't belong to certain irreducible representa- 
tion of the symmetry group for the undistorted cluster. This situation appears 
to be quite similar to the case of a single center, where the single electron dis- 
placement reduces the symmetry of system to the minimal one (axial). 

Supposing that other things being equal the configuration minimizing the 
inter-electron interaction has the lowest energy, we can introduce the following 
two-particle wave function 

^ro,o;r,(n, rl; q^^) = N (l ± A2) T«f i'J.,^ro7o(^"'i, 0) ^ro,o{r2, 0) , (25) 



where is the normalization factor, P12 is the electron permutation operator, 
is the operator of the symmetrized displacement gp- which transforms 

accordingly to irreducible representation r7 and operats in the space of the ith 
electron coordinates. The upper sign corresponds to the singlet wave function, 
the lower sign to the triplet one. The transformational properties r7 of the 
two-particle wave function (25) are defined by r7 = ro7o x [^7]^ for the singlet, 
and r7 = ro7o x r7 for the triplet. 



7 Electronic Jahn- Teller effect 

Non-rigid shell model gives a simple and obvious example of a local pairing 
within the two-electron ns^-like configurations to be a result of the correla- 
tion effects. The local pairing is promoted by the presence of a strongly polar- 
ized shell, as well as the orbital degeneracy or quasi-degeneracy within valent 
states (for simplicity, d-states) through the electric multipole s-d-interaction 
described by the effective "vibronic-like" Hamiltonian 

7 

where the V'' operator works within d-manifold, the are "vibronic-like" 
parameters. This interaction can result in a purely electronic Jahn- Teller effect. 

In general, one has to take account of the atomic displacements Q'* modes and 
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their interaction with electronic shifts: 

7 

This results in a complicated multi-mode Jahn- Teller effect with a correlational 
hybridization at the s-, d-electron modes and the local structural modes. 

This system will have all anomalous properties of a Jahn- Teller center, in 
particular, large values of the low-frequency polarizability. It appears that 
within the non-rigid shell model the completely filled electron shells do not 
quenched and can reveal many peculiarities similar to the nonfilled shells. 
The magnitudes of the shell shifts are correlation parameters which may be 
found by minimizing the energy functional E{q). The quantity A = £"(0) — 
E^qo) determines the pairing energy, i.e., the local boson binding energy. 

Non-rigid shell model can be considered to be a generalization of the well 
known shell model of the lattice dynamics and of the non-rigid anionic back- 
ground model by J.E.Hirsch ct al. [9]. In particular, a correlational pscudospin 
formalism can be successfully applied for a description of the valent states for 
the atomic systems with a correlational near degeneracy. Finally, we would 
like to conjecture the possible importance of the non-rigid shells correlation 
effects for a local pairing in copper oxides. 



8 Non-rigid shell model and hyperfine coupling 

A correlational shift of the one-electron shells could result in a considerable 
renormalization of the hyperfine coupling both for the nominally ns- and non- 
s-orbitals. Here we consider only two effects: 

1) an appearance of the effective non-s-contribution to contact hyperfine cou- 
pling, and 

2) an appearance of the effective nuclear quadrupole interactions for the nom- 
inally ns-electrons. 

Firstly, a shift of the one-electron non-s-shells for two-electron configuration 
implies an emergence of the effective electron density on the nucleus and could 
be detected in nuclear resonance, first of all in anomalous isotope (chemical) 
shift. In other words, the bare p-, d-, ... electrons within the non-rigid shell 
configurations appear to be involved in the contact hyperfine coupling. 

Secondly, a shift of the one-electron shells for two-electron configuration re- 
sults in a modification of the nuclear quadrupole coupling due to change in 
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the electric field gradient. Moreover, breaking the spherical symmetry of the 
electron density distribution within the ns^ configurations with the shifted 
shells leads to an appearance of the electric field gradient on the nucleus: 

Vij oc qiQj - ^q^Sij. 

The effects under consideration could provide real opportunities for a detection 
of the correlational shift of the one-electron shells with the help of various 
nuclear methods. 



9 Conclusion 

In a framework of the non-rigid shell model a correlated state of two-electron 
molecular configuration is described by a set of symmetrized shell shifts q'^ 
similarly to the well known shell model developed for a description of the 
lattice dynamics. Such a state could appear even for the completely filled 
shells thus resulting in a non-rigid atomic background with internal degrees of 
freedom. Contrary to conventional approach this background in common has 
nonzero electric and magnetic multipole moments. Non-rigid shell/background 
model results in a number of the novel unconventional effects including: i) 
a correlational mechanism of the local pairing; ii) a correlational (pseudo) 
Jahn- Teller effect provided by a joint account of the electron shell shifts and 
conventional nuclear displacements; iii) an appearance of the correlational 
current states. The model allows an introduction of the pseudo-spin formalism 
and effective "spin-Hamiltonian" for a description of the short- and long-range 
ordering of the non-rigid atomic backgrounds in crystals. Finally, the model 
can be readily built in the conventional band schemes. 
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Appendix 



Here an expressions of the matrix elements with a single-particle state 
il^k (r) = Nz,k e-^'- , where Nz,k 



^47r(2A; + 2)! 
are presented. 

The matrix element of kinetic energy of electron with functions of the same 
center is: 

t(a, a) = jirM f) (-|) U?) = (26) 



The matrix element of interaction of electron with the potential center with 
functions of the same center is: 

(27) 

I \ {2af 

2k + 2) n \ ' 

The overlap integral for the one-particle orbitals is: 

S[a,f3) = j dfipk{r- a) ipk{r- 
where p — Z\a — $\, 



u 



1 

a 



,-26 2k+l 



a 



El- 



1=0 



where a = Za. 



-p 2fc+2 



2j-l)\ I y ' \h) h\{a-h)\ ' 

- integer of s/2 , 



(=max{0,2j— n— 1} 

on+n'+a ^ ^ n+n'+2 

defined by — / (if|f - aPe-l^-"l|f - /^r'e-l^-^l = e"" J] p^C^"'"') . 

s=0 
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The matrix element of kinetic energy of electron with functions of different 
centers is: 



'2k+2 
I s=0 



2{2k + 2)\ 



2k+l 

E 

s=0 



2k 



+ 4A;(A; + l)5]p^Cf'*^-2) 

s=0 



(29) 



The matrix element of Coulomb part of inter-electron interaction is: 



iplin - a) ipl(f2-(3 



(if 1 (if 2 

fl - f2 

p {2k + 2)p^/ l\ 
2^k+\k + l){2k + 2)\ ^ ' 



(30) 



where 



G 



(fe) ^ 2' "^Yl (2fe + 2 - j) ^(2fcj-i) 

j=max{0,i-2fe-l} 



2^j\ 



The matrix element of exchange part of inter-electron interaction is: 



a(a,P^ = J Yp-y-^iJkin - a) ^k(f2 - 0) ^/'fe(fi - i)k{r2 - d) (31) 

4fc+5 k+l 2k+2 

where 



^i?(p) = -T E B\ 



— 1 
''s.n 



,■-"+1 



)(m,6-) (2j)! _ „ 
p2i+l ' 



min{j,m— s} 



S5(^) = sS(x)-^E,,„ (x,l); 



m— 1 



(2i + n + 1)! ~, 



2i+n+l 



s„„(x,e) = E ^^^":;r^^^ A^'"'^^ e 



=max{Oj— s} 



i=0 



X 



t=0 



t\ 



14 



p—x m—s 

2s+2j+ri (Y\ r l-l (cy \h 



In 27a; - ^0,/ - (-l)^+V-^z(-2a:) + (-1)^+" ^ ^n^-^M 



.=0 ^! 



min{i,m— s} s (O 

i=max{0,i-s+l} r=max{0,i4ti} V^' 'A 

tii^ ' l^/y'(m-/-s)!(2m-2/ + 2s + 2)! 

- the coefficients in '^af^'^^x^ — / (a;^ — t^) P2s{t)dt ; 
z=o _i 

the coefficients 6^*^ and D^^^ define the Legendre polynomial P2s(x): 
. ^ ^As) 21 As) (-1)'"' (2s + 20! 



2^^ (s - ly.is + 1)1(21)1 ' 
and the Legendre polynomial of second type Q2six)'- 

1 , X + 1 „ , , ^ ^(s)^2«+l . 



Q,^{x) = -ln--jP2six) -J2D\ 



1=0 

\l s-l 



rjis) ^ (-1)' V (-ir(4t + 3)(2t + 2/ + 2)! _ 

' (2/ + 1)! ^ 22*+i(2s -2t-l){s + t+l){t-l)\it + l + 1)1 ' 

7 = 1.78107 ; Ei{x) - the exponential integral function [1]; 

^ a + (3 a — P 

C = — ^ , ^ 



\a — (5\ \a — (5\ 



The matrix element of interaction of electron with the potential center with 
functions of different centers is: 



u 



(dj) - j -li;,{f-d) ^fe(f-/3) (32) 

„2fe+2 k+1 ( 2k+2 

^ ^(2^T2)! + ms{ri)\Q2s{v) E - e-^^lpOi ^if/'Hp) 



g-P? ~ (fc+i) 



P 

2s+2i 



+ 



s+2i ^KS) /I \ , 

E ((-l)^+^e''£;z(-p(e + 1)) + e-^Ei{-p{i - 1)) 



i=0 



Z=0 P 

'-1 ^/i 



^-"^ E ySn,i {{-i)'^\i + 1)' + (e - 1)'^))] } ■ 
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Tables 



Table 1. The results of numerical minimization of the full energy functional 
for the Is^-configuration {k — 0) a,t Z — Zq. 



Zj 






f 




1 

i.O 


n 078 
U.U I o 


n 078 
U.U I O 


O. i4iO 


— i.Oi4U 


i.O 


U.UOo 


U.UOo 


O. i4iO 


— i.00i4 


1 7 


u.uoy 


u.uoy 


1 /1 1 

O.i4i0 


1 8987 
— i.oZo / 


i.O 


U.UOi 


U.UOi 


o.i4iO 


— Z.iiD 


1 Q 

i.y 


U.U4D 


U.U4D 


o.i4iO 


— Z.4ZOD 


2.0 


0.043 


0.043 


3.1415 


-2.7510 


2.1 


0.041 


0.041 


3.1415 


-3.0984 


2.2 


0.038 


0.038 


3.1415 


-3.4658 


2.3 


0.036 


0.036 


3.1415 


-3.8533 


2.4 


0.034 


0.034 


3.1415 


-4.2608 


2.5 


0.033 


0.033 


3.1415 


-4.6882 


Table 2. 


The results of numerical minimization 


for the 1 


s^-configuration [k = 0) at Zq = 2.0. 


Zj 


a 




f 


E{dJ) 


i.O 


0.0 


0.0 




-2.8125 


i.D 


0.0 


0.0 




-2.8400 


1 7 
i. / 


0.0 


0.0 




-2.8475 


1 8 

i.O 


0.0 


0.0 




-2.8350 


1 Q 

i.y 


0.011 


0.011 


3.1415 


-2.8026 


2.0 


0.043 


0.043 


3.1415 


-2.7510 


2.1 


0.065 


0.065 


3.1415 


-2.6809 


2.2 


0.084 


0.084 


3.1415 


-2.5931 


2.3 


0.102 


0.102 


3.1415 


-2.4877 


2.4 


0.118 


0.118 


3.1415 


-2.3650 


2.5 


0.132 


0.132 


3.1415 


-2.2253 



of the full energy functional 
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Table 3. 



g-expansion of the matric elements up to q^. 



General expression 



Q'-expansion up to g 



S {q, 



1 - g2 2 / dfP 



1 



,2 2Z^ 



3(2fc+l) 



Z2 



2(2A:+1) 



u{q,q)^jfa'' + q^jf{b'' + ac) 



Z-q^^, k^O 



t {q, —Q) ~ — H-/^ ^^"^ Aa + q'^ J dfc A a] 



_ 2 

2(2fc+l) ^ 4fc^-l ' 



A; 7^ 



u {q, -q) « / fa' + q' I f {-b' + a c) 



Z^q'^, k = Q 



z _q2 



2Z3 



fe+1 



3(2fc+l)(fc+l) ' 



c -q) ~ / ^alal+ 



Z(^- 

'12 ^ ^ 

+g2 2 / (-201610262 + afbl + a?a2C2) -g'Z 



(4fc+3)! 
Vfc+1 24^+2 [(2fc+2)!]^ 
2y3 (4fc+2)! 



3-2-*fc-i[(2fc+2)!1^ 



a{q,-q}^I^alal+ 

+q'2J^{-albl + ala2C2) 



y (4fc+3)! \ , 

\k+l 2«+2[(2A:+2)!lV ~^ 



2«+2[(2A:+2)!] 



3(2fc+l)(fc+l) 3-24fc-i[(2fc+2)!]2, 



(4fc+2)! 



£;(o) = - / dfaAa - 2Zo f ^a'+ 



d£idf2„2 „2 
l"2 



+ 1 "^t^a^d' 



Z'-2ZoZ + lZ. 

Z2 _ 2ZZo 



k^O 



2k+l fc+1 



+ 



(4fc+3)! 



7 ^ 

^ Vfc+1 24fc+2 [(2/0+2)!]'' 



A; 7^ 



£;(2) = - J dfcAa- J dfaAaJ dfb'^- 
-2Zo{jfac + jfa'Jdrb') + 
+2 / ^ [a?a2C2 - 01610262 + afai / dfb'] 



4Zf_ 
3 



^4 4(fc+l) 
^ 3(4fc2-l) 



Z — + ^ j , k 



-Z3 



(4A;+3)! 



3-24'=+i(2fc+l)[(2fc+2)!]^ 



ky^O 



Here ri2 = I ri — r2 I, and the index i of the functions a, 6, c indicates depen- 
dence on fj. 
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Table 4. 



k 


E ■ 




^ . 




^7 f r\\ 7 





-2.84766 


0.0 


1.6875 








1 


-2.22965 


0.3437 


2.7110 


0.2205 


-0.2559 


2 


-1.79140 


0.5061 


3.2084 


0.3463 


-0.5204 


3 


-1.48131 


0.6385 


3.4999 


0.3638 


-0.7030 


4 


-1.25036 


0.7644 


3.6778 


0.3424 


-0.8192 
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Figures 




Fig. 1. The displacement vectors a and /3 define the centers of the one-particle 
orbitals relatively to the potential center. 




0.0 0.5 1.0 1.5 



2.0 



2.5 3.0 



r, a.u. 

Fig. 2. The Slater orbitals ^/;(r) (6) at Z = 



(fc) 

rain ' 




Fig. 3. The full energy functional E{q, Z) at Z = Z^j^„. 
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O.O' «' -'-T-T I 

1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



Z 

Fig. 5. The value of the electronic shell displacement q minimizing the full energy 
as function of the given parameter Z. The points correspond to minimal values of 
the full energy at the given k. 
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0.04 H 




■O ^^0 



Fig. 6. The two-electron density distribution in plane 2 = for the state 
^ (r*!, 9; — A; = 1, = 2 a) without shifts of the electronic shells (Z = Z^^j^ 
at g = ); b) with shifts Z = q = (g|^L,0,0). 
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